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1. Three-Nucleon Interaction V tr " . 

In the description of the Urbana potential series [2], the three-nucleon interaction is presented as a sum 
of a long-range part derivable from the two-pion exchange diagrams, and a phenomenological short-range 
part, adjusted to reproduce the binding energy of the three-body nuclear system. The two-pion exchange 
interaction (V 27Tl tni) is given as 

V^tni = A 2ATv,T 2 ,T 1 -T 3 }{S 12 T(r 1 2) + cTi-<T2Y(r 1 2)),S 13 T(r 13 )+a 1 -a 3 Y(r 13 ))} 

cycl. (1.1) 

C 2 An-,T 2 ,T 1 ■ T 3 ][S 12 T(r 12 ) + ^ ■ a 2 Y(r 12 )), S 13 T(r 13 ) + ^ ■ a 3 Y(r 13 ))} 

Here X) C ycZ represents a sum over cyclic permutations over the indices 1,2, and 3. r, a, and SV, are the 
isospin, spin, and tensor operators, and {,} and [,] denote the anticommutators and commutators. The 
T(r) and Y(r) are radial functions associated with the tensor and Yukawa parts of the one-pion-exchange 
interaction, and C 2lx = \ A 2lx . 

The short range repulsion is phenomenological and is given as 

Vfl, tm -[/o^T 2 (r 12 )T 2 (r 13 ) (1.2) 

cycl. 

Note. For the Urbana IX potential, the fitted parameters are: A 21T = —0.0293 and Uo = +0.0048. 

The relevant matrix elements for our calculation are of two types only. The matrix elements of the 
form V h ™ ,a „ .. . h are derived in section 2, whereas the matrix elements of the form —V h ™ ,a „ . h h are 
presented in section 3. 

2. Density-Dependent Matrix Elements. 

In this section we work out the computation of the integrals 

Vta\, bl ,P,a 2 M =< ^(1)^(2)^(3)1^1^(1)^(2)^(3) > . (2.1) 

We employ the same methods as developed for the calculations of the two-body matrix elements, namely 
using Fourier-transforms in order to separate the variables. We use an expansion into Harmonic oscillator 
functions H^(r) for which the Fourier transforms are again Harmonic oscillator functions H^(q). 

We will make use of the particular angular momentum coupling in order to make the computations 
feasible. For the matrix elements given by (2.1) and requiring jh = j p the angular momentum coupling is 
similar to that of our two-body matrix elements: 

< (aia 2 )A|U e ^|(6 2 6 1 ) A >=S mh ,m p (-)*•* +fc * (-)(ka 2 +k bl+ m a2 -m bl ) 

< j ai m ai j a2 - m a2 \\n >< jb 2 m b2 jb 1 - m bl \\fi > V mh 
< (aia 2 )A|U e ^|(6 2 6 1 ) A >=6 mh , mp (_)*.!+** (_)(*.* +*»i+™. 2 -n» tl ) 

< i ai m ai j a2 - m a2 \\n >< 

jb 2 TTT'b 2 jb 1 W^jjA/i > ( )V mhmaiTnbi ,m a2 m p mi, 2 

(2.2) 

which includes the "Ring" -phase for two-body p/i-matrix elements. Here we do sum over — m p and 
we do sum over all other m's. This angular momentum coupling applies as well for the matrix elements of 
section 3. In this section we handle the sum over the cyclic permutations by calculating the three matrix 
elements separately 

V^ uhmM = < ^(l)0 ai (2)0 fcl (3)|y 4 -|0 p (l)0 a2 (2)0 fc2 (3) > 

< ^(2)^(3)^(1)1^1^(2)^(3)^(1) > (2.3) 

< ^(3)0 ai (l)0 fcl (2)|\/ 4 -|0 p (3)0 a2 (l)0 fc2 (2) > 
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2a. Short range repulsion term. 

As this term does not contain any additional operators it is particularly simple to handle. We use Eq. (2.2) 
of Reference [3] 

T 2 (r l2 ) = 4tt - f q 2 dqf(q)^(-yij e (qr 1 )je(qr 2 )[Y i (h)®Y e (r 2 )f ) (2.4) 

T(q) = J r 2 12 dr 12 T 2 (r 12 )j Q (qr 12 ) (2.5) 



using 



Correspondingly we write 

T 2 (r 12 )T 2 (r 13 ) = (4tt) 2 £ (-)^H 2 ! 3 

t Sdq 2 f(q 2 )je 2 (q 2 r 1 )je 2 (q 2 r 2 ) x J qtdq 3 f{q 3 )j l3 {q 3 r l )ji 3 {q 3 r 3 ) ^'^ 

x [Y e2 (h) K fa (f 2 )] (0) [y, 3 (n) ® ^s(r 3 )] (0) 
recoupling the spherical harmonics results in 

T 2 (r 12 )T 2 (r 13 ) =(4t) 2 4= ]T 



47T 

<?l/ 2 /3 



X 



^ y qldq 2 f{q 2 )j l2 {q 2 r 1 )j l2 {q 2 r 2 ) X ^- J qldq 3 f{q 3 )j l3 (q 3 ri)ji 3 {q 3 r 3 ) ( 2 - 7 ) 



re m (°) 

y £l (f!)® [y £2 (f 2 )®r, 3 (f 3 )] (£l) 

This leads to l\ = for the first matrix element, l 2 = for the second and ^3 = for the third. In turn we 
can write the interaction for the three matrix elements respectively as 



T 2 {r l2 )T 2 {r l3 ) = ^R £ i [V^Y t (f 2 ) V^Y t (f 3 )] 

£ 

x- [qldq 2 f(q 2 )j e (q 2 n)je(q 2 r 2 )x - I ' qldq 3 f(q 3 )j e (q 3 n)je(q 3 r 3 ) 

7T J 7T J 



(2.8) 



for the first matrix element and 

T 2 (r 12 )T 2 (r 13 ) = £ [v&/(fi) V^T/fa)] 



x ~ J ( &dq2f{q 2 )j a {q 2 r 1 )j (q 2 r 2 ) x ^ J lld^iq^jfXq^jfXq^) 
T 2 (r 12 )T 2 (r 13 ) = ]T [V^Y e (h) V^Y e (f 2 )] 

£ 

x ~ J <l3 d< tiT(q 3 )j a (q 3 ri)j (q 3 r 3 ) x ^ J qldq2T(q 2 )ji{q 2 ri)ji>{q 2 r 2 ) 



(2.9) 



for the second and third respectively. 

For the first matrix clement we write 



< 



^(l)0 ai (2)0 fcl (3)|T 2 (r 12 )T 2 (r 13 )|^(l)^ 2 (2)0 b2 (3) >= R A A J rUnRh^R^n) 

x ^ <i ai ||y A ||ia 2 > J r 2 2 dr 2 R ai {r 2 )R a2 {r 2 )^- J qldq 2 f(q 2 )j x (q 2 r 1 )j x {q 2 r 2 ) 



x(-) fe <>2^ < jb 2 \\Yx\\j bl > J rjdr 3 R bl {r 3 )R b2 {r 3 )?- J q 2 3 dq 3 f(q 3 )j x (q 3 n)j x (q 3 r 3 ) 

(2.10) 



We now assume R ai (r)R a2 ( r ) is expanded into harmonic oscillator functions as 

R ai (r 2 )R a2 (r2)=J2 A n H n(r2) (2.11) 

n 

Using Gauss integration the coefficients can be written in terms of the functions at the Gauss-points, r i; 
with weights, Wi, as 

K =J2 R MRa 2 (r i )H^(r i )^ Wi (2.12) 

i 

where the barred functions are those without the exponentials. In turn, we can write the integral 

/ rldr 2 R ai {r 2 )R a2 (r 2 )- [ qZdq 2 f(q 2 )j x (q2ri)j\(q 2 r 2 ) = 

J n J (2 13) 

where 

T*, n = f q 2 2 d q2 H^q 2 )f(q 2 )H^ q2 ) = £ q 2 w s H^q s )f(q s )H^q s ) (2.14) 

We proceed similarly with the integration over #3 and r 3 . Thus this matrix element can be written as 
< ^(l)</» ai (2)^ 1 (3)|T 2 (r 12 )T 2 (r 13 )|^(l)<A„ 2 (2)<A 62 (3) >= 

^<j ai m\\ ja2 >(-) k ^<. 

A A 
^2 Ra t (n)R a2 {ri)Rb t (r 3 )Rb 2 (rj)AS l 



{-) kb ^ <j ai \\Yx\\ja 2 > (-) kb ^<j b2 \\Yx\\j bl > (215) 



X 



where 

AS itj = H^rjwMir^w^T^D^ (2.16) 

with 

D X m ,i = (-) A A / r 2 dr s H^(r s )H?(r s )R p (r s )R h (r s ) = (-) x \Y,r 2 w s H^(r s )H?(r s )R p (r s )R h (r s ) (2.17) 

■> s 

For the second matrix element we obtain 
< ( /) 6l (l)^(2) ( / )Ql (3)|T 2 (r 12 )T 2 (r 1 3)| ( /» 62 (l)0 p (2)0 a2 (3) >= J r 2 2 dr 2 R h (r 2 )R p (r 2 ) 



x(-)*° l -jp < 3bA\Y\\\j b2 > J rldnRb^r^R^in)^ / <iidq 2 T(q 2 )./oU/2>-j. )Mw-2 ) 



A-? a *^Y < 3a 2 \\Y X \\j ai > J r 2 dr 3 R ai (r 3 )R a2 (r 3 )l J qldq 3 f{q 3 )j x {q 3 n)j x {q 3 r 3 ) 

(2.18) 

Again, this equivalent two-body matrix element takes the form 

< fcl (l)^(2)0 ai (3)|r 2 (r 12 )T 2 (r 13 )|^ 2 (l)^(2)0 a2 (3) >= 

V4tt . . ... k \/47r 

^<3bA\Yx\\3b 2 > (-) ka2 ^ <; 
A A 

K (n)Rb 2 (n)R ai {r )Ra 2 (rj)AS^ 



(-) kb ^<3b 1 \\Yx\\jb 2 > R fe -^ <3a 2 \\Yx\\j ai > (219) 



In this case 



where 



C„° - / r 2 2 dr 2 R h (r 2 )R p (r 2 )H° n (r 2 ) = Y / r 2 k w k R h (r k )R p (r k )H° n (r k ) 



Similarly, the third density term is 

T 2 (r 12 )T 2 (r 13 ) = Q 3 E [Yt(h)QY t (f 2 )] 

I 

X Jq 2 dq 2 f(q 2 )j e (q 2 ri)j e (q 2 r 2 )x j qldq 3 f{q z )j () {q z r 1 )j {q 3 r 3 ) 

The form can be obtained from the previous case by using the transpose i.e. ASf j — AS 1 ^. 

2b. Two-pion exchange term. 

For the model V 27r-exchange TNI takes the form 

y 2 *™>° : =:A2 7r {Tir2,Tir3}{(5 , i2r(ri2) + aia 2 Y(n 2 )), {S 13 T(ns) + <7i<7zY(r 13 ))} 
We can write the tensor interaction as 



r W (r 12 )® [(7i®<7 2 ] 



(0) 



S 12 = 3 < 10101X0 > V4tt 
with K = 2, whereas the sigma interaction is identical with K = 0. By defining the form factors 

4tt V K (q) = An J V K {r)j K (qr)r 2 dr 
with V(r) = y for if = and V(r) = T for K = 2 the interaction can be written as 

A 2n (4iT) 2 ^2Y^ < 1010|i^ 2 >< 1010|X 3 > i- < 40^0|X 2 >< 4040|if 3 > 
x 2 if 3 X2 ^ 3 



1 1 g 2 2 ^2^fe)iii 2 (*) ( ' 1 - £2 - K2) ^ 1 fer 1 )^ 2 ( < Z2r2) 
5 3 2 ^3^ 3 (93)44(i) ( ' 4 -' 3 - X3) i£ 4 (?3ri)j£ 3 (?3r3) 



[^(fl)®!^)]'* ® [<7l®<7 2 ] 



l(^2) 



(0) 



[y 4 (f 1 )®y 3 (f 3 )] (if3) ® [a!®a 3 



i(*s) 



1(0) 



(2.20) 
(2.21) 



(2.22) 



(2.23) 



(2.24) 



(2.25) 



(2.26) 



Here K 2 or K 3 are either 0, for the a ■ a-term or 2, for the tensor term. These need to be recoupled: 



[Y^fi)®!^)] yK2) ® [<71 ® <7 2 ] 



1(^2 



(o) r 



[y 4 (f 1 )®y 3 (f 3 )] (K3) ® [ffl®<7 3 ] 



1 (^3) 



(0) 



(2.27) 



[Yt^h) ® (Ti] (Xa) ® [Yt a (r 2 )®<r 2 



(A 2 



(0) 



[Y e4 (h)®<Ji] (X3 ' ® [Y ea (r 3 )®<T3 



(A 3 



(o) 



This, in turn we can write as 



J 

2 A3 



=(-)«-(-)«-A 2 A iJ f S A3 5 4-{' 1 ' t t ID 

1 (J) r 

[y 2 (fa) ® a 2 ] (A2) ® [y 3 (4) ® a 3 ] (As) 



[^(ri)®^]^® [y^(fi)®<7i] 
£1 1 A; 



-i(o) 



= {-Y^ +l ^k 2 k 3 x 2 x 3 m^l 1 A'jj^ 1 ^ 2 f 2 }{ £ ! l { 



[Y tl (fi) ® ^4(^1 )] L O [cti ® cti 



(S) 



(J) 



[y 2 (f 2 )®a 2 ] (Aa) ® [y/,(f 3 ) ®<7 3 



(A3) 



i (o) 



(2.28) 



= (-) 



K 2 K 3 X 2 X 3 JhhS 



tx 1 A 2 

^4 1 A3 

L S J 

(J) 



J I 1 1 A 2 J I 1 1 A3 J 



[y 2 (f 2 )®<7 2 ] (A2) ® [y 3 (f 3 )®a 3 ] 



(a 3 ) 



t(0) 



(2.29) 

As the reduced one-body matrix element of (Yg ® cr) A vanishes for A = 0, only the first matrix element in 
(2.2) gives a contribution. For it we have 5 = and L = 0. This implies J = and A 2 = A3 as well as 
l\ = £4. This term leads to a density dependent finite range Migdal g or tensor interaction. For this case 
we evaluate the 9-j symbol and simplify the interaction to 



£1 
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£2 K 2 
1 A 



h K 3 



} [[y 2 (r 2 ) ® <r 2 ] (A) [y 3 (f 3 ) ® <r 3 ] (A) ] ^ (2.30) 



By using Eq. (2.1) of Reference [3] we write the matrix element as 



< (aia 2 )A|V e// |(6 2 6i) A >= 9A 27r 

53 53 ^2^3 < 1010|X 2 >< 1010|X 3 >< £iO£ 2 0\K 2 >< £iO£ 3 0\K 3 > 



K 2 =0,2 K 3 =0,2 



^3) ( _ )( , 1+ A +1 )^{4 £2 Jfcjj*! t 3 X 3 | 



/ 

2 

7T 



r?driiJ h (n)i2p(ri)^(-) fc -i < j ai \\[Y^v) x \\j a2 > ±-^(-)^ < j„ 2 \\[Y^a] x \\j bl > 



qldq 2 V K2 {q 2 ) J r\dr 2 j^ {q 2 r\)u 2 (q 2 r 2 )R ai (r 2 )R a2 (r 2 ) 
^ j qldq 3 V K3 (q 3 ) J r^r 3 ^ 1 (g 3 ri)^ 3 (g'3r 3 )i?6 1 (r 3 )i? b2 (r 3 ) 



(2.31) 

Exchanging variables 2 and 3 results in the identical expression. Thus, the commutator vanishes for this 
density dependent term, whereas the anti-commutator obtains a factor of 2. Further, we use {t\t 2 ,t\t 3 } = 



2-T2T3. Again, we write the integrals as sums over Gauss-points 
< (aia 2 )A|^ e// |(6 2 6 1 ) A >= 4 • 9A 27r < f 2 f 3 > 

t&Ma < 1010|if 2 >< 1010|if 3 >< £iO£ 2 0|ir 2 >< hOl 3 Q\K 3 > 



K 2 =0,2 K 3 =0,2 

K 2 \(h l 3 K 3 



{j){ t 2 -K 2+ t 3 -K 3 ) { _ ){ X + l)^l *2 /C2j|«l h A.3j (2 32) 



4?r (-) fe - <ia 1 ||[^ ) ^] A |b 02 > ^PH*" <^II[^ 3) -] A |I.^ > 



&rk,,V£V tlt >V* M H£ (n)rf Wl Ht (r 3 )r*w 3 R ai ( n )R a2 (r t )R bl (r,)R b2 fa) 
where we have used the same definition of D as before (3.17) and 

V™ =^2ll^S e MS e n 2 (qk)V K (q k ) (2.33) 

The selection rules require l\ + £ 2 = even and l\ + £ 3 — even which in turn requires £2 + ^3 = even. Further, 
£1,^2, £3 are all restricted to values of A, A ± 1. 

We now discuss the four cases separately. The first case, K 2 = 0,K 3 = leads to a density dependent 
(7i <72 interaction: 

< (aia 2 )A|V^|(Mi)A >= 4A 27r < t 2 t 3 > D'V^fHiir^WiH^y.Wi 



<j ai \\[Y^a^\\ ja2 > ^{-) k »> <j b2 \\[YWaf\\j bl > (2-34) 
A A 

R ai (n)R a2 {n)R bl {rj)R b2 fa) 
second case (i^ 2 = 0,K 3 = 2): 

< (aia2)A|^ tensor |(62&i)A >= 4A 2 , < f 2 f 3 > >/6M 3 (i) (/s+/l) < *i040|20 > { \ \ \ } 

— — (2.35) 

47r (-) fe - < jaJ[^V] A ||j a2 > ^(-)^ < ife2 ||[F^) (T] A|| i6i > 



A A 
# ai )Bo2 fa)-Rf>i fa)-Rb 2 fa) 

third case (if 2 = 2,^ 3 = 0): 

< (aia2)A|^ e/ ^™|(626i)A >=4A 2 . < f 2 f 3 > V&hh < *iM 3 0|20 > « (/s+/l) { ^ t J} 

^ iS y^ 1K o,«3^ 1(ri)r 2 Wi ^3 (r . )r 2 w . 

, — — (2.3b) 

477 (-)*•! < jaJ[^V] A ||j a2 > ^(-) fc * < j b2 \\[Y^a] x \\j bl > 



A A 
fioj {n)Ra 2 {n)R ai {r 3 )Ra 2 fa) 

This term is similar to the previous term. They represent a density dependent tensor interaction. The last 
term, the tensor squared term has (K 2 = 2,K 3 = 2): 

24A 27r T 2 f 3 hhhh < 40^ 2 0|20 >< £i040|20 > 
W (^3) ( _ )( A + i){^ to 2J|4 /, 2J 

(2.37) 



[n H ka > < j ai || [Y^a} x || j a2 > ^(-)^ <.7 fcl ||[r^) ( 7] A || Jb2 > 



A A 

^, s ^',l l£2 K 2 / l£3 ^ 2 (n)rf Wi Ht fa)r>,- i? ai fa)i? a2 fa)i? 6l fa) j? &2 fa) 
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From this case we can split up terms that are similar to the previous ones. We write 



24A 27r < f 2 f 3 > 4444 < 4040|20 >< 4040|20 > 



(2.38) 



-4^(-) fe - <j ai \\[Y^a] x \\j a2 > ^L{-) k ^ <j b2 \\[Y^a] x \\j bl > 
A A 

#ai (n)-Ra 2 (n)-Rfei {rj)Rb 2 ( r j) 

The remaining term is: 

< (aia 2 )A|l /e// ' fc = 1 |(6 2 6i)A >= 72A 2ff < f 2 r 3 > 



-S(-)'- <i»ll[v"V] A iii«>^(-)'>. <i tl ii[y" l »] A iiA, > 

A A 

DtsV^Xf^H^yi^tir^wj R ai (n)R a2 (n)R bl (rj)R b2 (rj) 



(2.40) 



-^H*" 1 <ia 1 ||[^ 2 V] A ||i a2 > < j b2 \\[Y^a] x \\j bl > 

A A 

^^^^^^(rOrJttfi^^Or^j R ai (n)R~a 2 (n)R bl (r 3 )R b2 (r 3 ) 

The term with fc = can be combined with the term (2.34). Also, the term with k — 2 can be combined 
with the terms (2.35) and (2.36). 

Thus, we end up with three terms. The first corresponds to a density dependent o\<Ji term and can be 
added directly to that term with the form: 

< (aia 2 )A|^ e//> |(6 2 6 1 ) A >= ^ < t 2 t 3 > HfoitfwiHfcjtfwj 

il 

( _ )( m + i)^R feoi < j ai \\[Y^\\j a2 > ^(-)^ < jb2 \\[Y^an jbl > 
A A 

R ai (n)R a2 {n)R bl {rj)R b2 (r,) 

The second corresponds to a density dependent tensor interaction: 

< (ai«2)A|^ e/Xtensor |(&26i)A >= 4A 2n < t 2 t 3 > H^in^WiH^ir^rjwj 

D^V^Xf 3 " 3 + D% S V^V:;!^ + 5^(21 + !)(-)* 

e 

< 1040|20 X 1040|20 > f4 4 21 | 2 2 2 1 £ V 2,« 1V 2,/ 

< 4040|20 > 12 2 £/ll 1 1J ™,s V m,n V s,t 

Veiih < 4040120 >(i)(^o{^ 2 1 



(2.41) 



(2.41) 
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3. Exchange matrix elements. 

In this section we work out the exchange matrix clement which we define as 

ytni,x ytni.a 

pihi;h,2P2 h,pi,p2\hi,p,h2 

We assume the interaction can be written as sum over terms each having the form: 

-T^\l)®[T^®T^] (Xl) f ] 
This leads to the exchange matrix element as 

KT.MP2 = "{ < Ml)Pi(2) P2 (3)|[r i (Al) (l)® [7f 2) ®7f 3) ] (Al) j ( Vi(l)p(2)^(3) > 



(3.1) 



(3.2) 



+ < h(2) Pl (3)p 2 (l)\ if ® [if 2) ® T. 



-,(A 3 )-|(Ai) 



(0) 



|/ ll (2)p(3)/ i2 (l) > (3-3) 



+ < M3>i(l)P2(2)|[7f lJ (l)® [T 2 (A2) ®T 3 (A3) ] (Al) J^>i(3Ml)^(2) >} 

Here the three terms arise from the cyclic permutations. Using our phase convention for ph-ph angular 
momentum coupling given by Eqs. (1.4, 1.14) of Reference [3], we can carry out the summation over all 
m's except m p = mu, using the Wigner-Eckart theorem we obtain: 



(0), 



- r^+^+ a ^+ a >^{ Al A2 A } x {i<^ilfci>i<pill^-ii^>i< Pa ||23Mi^ 

A L Jpi Jhi jh > y Ai A2 A 



) 



A 



Ai A 2 

f (p 1 ||T3^||/ 1 )i(p 2 ||T A ||/ l2 )J-(/ l ||T A 1 ||/ ll )} 
A2 A Ai - 1 

(3.4) 



3a. Short range repulsion term. 

We first turn to the correction term due to the short range repulsion as the most simple contribution of 
ytm_ p rom sec tion 2 we take the interaction as 



T 2 (r 12 )T 2 (r 13 ) =(4tt) 2 -^ V t 2 £ 3 < * 2 M 3 0|*iO > 

x ^ y ( & dc tiT{q2)3i2 {l2ri)je 2 (<? 2 r 2 ) x ~ / lldq 3 f(q 3 )j e3 (q 3 ri)j/> 3 {q 3 r 3 ) ( 3 - 5 ) 

x "^(fi)® [^ 2 (f 2 )®y £ 3(f 3 )]^ l) ] (0) 

Combining this with Eq. (3.4) and considering the A's were renamed in order to obtain (3.4) we find: 

v tni,x,\ _, N(fc„ 1+ fc fa + fc,„+A) ( 2£ 1 + !)( 2 ^ 2 + 1 ) f h jh t jh \ tfn/nlxn 

•Wi,h2j* A I i P1 ^ 2 A J £lUt2U i AU 



( _)*, a xiE < P2 ||y A ||/ l2 >] [(-) fe ^ </*||y/j/n >] [(-) fe "^ <Pi||^ 2 b> 

A 



> 

47T 

2 



xjy rfdnRhir^Rh^n) 



q 2 dq2T(q2)j l2 {q 2 ri) / r 2 dr 2 R pi (r 2 )R p (r 2 )je 2 (g 2 r 2 ) 



^y 93 dq 3 f(q 3 )j\{q 3 r 1 ) J r 2 ,d3 2 R P2 (r 3 )R h2 (r 3 )j\{q 3 r 3 ) 
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^ J qldq 2 f{q 2 )j tl {q 2 n) J rldr 2 R hl (r 2 )R h {r 2 )j^{q 2 r 2 ) 
^ J qidq 3 f(q 3 )j e2 (q 3 r 1 ) J rld3 2 R Pl (r 3 )R p (r 3 )je 2 (g 3 r 3 ) 
+ J r\driR pl (ri)R p {ri) 

^ J qldq 2 f(q 2 )j\{q 2 n) j rldr 2 R P2 (r 2 )R h2 (r 2 )j\(q 2 r 2 ) 
\ J qldq3f(q 3 )j ei (q 3 r 1 ) J rld3 2 R hl (rs)R h (r 3 )j ei (q 3 r 3 )^ 
Using Gauss integration, we write the radial matrix element as 

(r k )T*R P2 ( rj )R h2 ( rj )H? {r )r 2 jWj 
+r 2 w k R p (r k )R Pl {r k )H% (r k )T^ n R hl (n)R h (n)H 1 ^ ( n )r 2 Wl 
(r k )T s \R P2 (rj )R h2 (r, )H* (r, )r 2 Wj 
+r 2 w k R P2 (r k )R h2 {r k )H% (r k )T^ n R hl {n)R h { n )H^ {n)rjw t 
H 1 / (r k )T^R Pl {r 3 )R p {r 3 )H^ ( rj )r 2 Wj } 

defining the kernel 



(3.6) 



(3.7) 



G l (r k ,n) = yJ^H^r^Hfcrfwi (3.8) 
allows us to write the radial integrals as 

[Rhir^Rh^KA^Mr^irk.r^ir^r,^ 

+fi P (r fe )ii Pl (r fe )^ 1 (rO^(rOG^Hr-fe,ri)G ,A (r fe ,r J 0fi P2 (r^ J R/ l2 (r- j ) (3.9) 
+R P2 (r k )Rh 2 (r k )R hl (n)R h (n)G^(r k , n)G l * {r k , rj )R Pl (r,)i? p (r,)} 
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3b. Two-pion exchange term. 

Again, we take the interaction from section 2 as 

A 27T (4tt) 3/2 < 1010 l^20 >< 1010|A 3 > (-J^+^s+As 

K 2 =Q,2 K 3 =0,2 

- I qldq 2 V K %q 2 )i 1 i 2 (tf 1 - e2 - K2 \ ]ll (q2ri)h 2 (q 2 r2)^- < £iO£ 2 0|A 2 > 

7T 

2 

7T 



K 2 

qldq^v K3 {qz)iMi) {li ^ K3) 3iAq^i)jiAq^)4r < 40£ 3 0|A 3 > 

A 3 



^i(ri) <g> [o-i <g>CTi 



(Ai) 



[r, 2 (r 2 )®a 2 ] (A2) ® [F, 3 (f 3 )®a 3 



(As) 



(AO- 



CO) 



(3-8) 

Here A 2 or AT 3 are either 0, for the c-cr-term or 2, for the tensor term. Further, 5 is in the anti-commutator 
term and 1 in the commutator term. As for the present correction we only have a contribution from the 
anticommutator, we write the interaction as twice the 5=0 contribution: 



18A 27r (47r) 3 / 2 < 1010 l^20 >< 1010|A 3 > (- 

K 2 =0,2 K 3 =0,2 

2 



.A2-M3+A1 



qUq2V K %q2)£ii2(if 1 - i2 - K2) jt 1 (q2ri)je 2 (q2r 2 ) < £i0^0|A 2 > 
\ J qtdq3V K3 (q 3 )iJ3(i) {e4 - e3 - K3) jeAQ3n)j i3 (q3r 3 ) < h0£ 3 0\ K 3 > 



(3.9) 



^MZ I IKI T till t 



[Y, 2 {f 2 )®a 2 ] {X2) ® [Ye 3 {r 3 )®a 3 



(As) 



(Ai) 



1(0) 



We carry out the summation over all m's and do the summation over cyclic permutations by using Eq. 
(3.4). However, it should be noted that in the second and third term of the cyclic permutations the A's in 
Eq. (3.4) were renamed 



18A 27r < 1010|A 2 >< 1010|A 3 > Ixldxlrfds 



AiA 2 r Ai A 2 A 

A Jpi 3h± Jh 



{ 



}(-)** 

^(ii-e 2 -K 2 +e,-e 3 -K 3 ) < i^q^q >< £ 4 040|A 3 > (_)( fe »+ fe P i+ fe >n) 
^ J q 2 dq2V K2 (q2)ji 1 (q2n)j£ 2 (q 2 r 2 ) 
^ J qldq 3 V Ka (q 3 )je 4 {q 3 ri)ji 3 (q 3 r 3 ) 



A 



[(_)feft ^ :</ ' l|:rAl|l ^ 1>] [f-)^ 1 ^ 1 ^ 1 !!^^^!^] [(-) fep2 ^x i: <^l|[^3^ A ll^>] 



A Ai A 2 
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[(-)*»^(y.r Al lb)][(-)*"^fell[v; 1 »] A l|A2)][<-)'^(''IIK.<'] A 'll''i>]} 

We write the radial integrals as summations over the Gauss points. We compute the matrix element as 

ty(ti-t a -K a +u-t a -K 3 ) < ^ ^ 2 o|if 2 >< e 4 0£ 3 0\K 3 > (-)^+ k Pi+W) 



[(-)'»^{/,||F»MlM[(-)'"^(wllK,"] A 'b)][(-)'"^(KllK,<'] A llM 

Ai A2 A 



[(-)*''^<Bi|y A iift 2 )][(-)* , '^<'>ii[5fe»] A MiAi>][(-)*"'^(piii|y,,<T] i -iip>] 

A Ai A2 

r^ fe i? p2 (r k )R h2 (r k )H^(r k )T^ ^ R hl (n)R h (n)H^ (r^w, 

£ iH? t till 1 t ^}(^ 4 o,A 2 o >( -)-3 



[ ( _)^ 1 ^I ( p 1 ||yA 2|b)] [ ( _ )fcp2 V^ (p2 || [y£2(T] A|| M][( _ )fe ^^||^ 3(T] A 1||fti)] 

A2 A Ai 

^^ 1 (r t )J? p (r fc )ffltOi' l4A ^ 1 (^iji.('-.)^(r,)^ 1 

Hi 1 (r k )T s f' e ^R P2 {r 3 )R h2 {r 3 )H^ ( rj )r 2 jWj 
We define the interaction kernel as: 

K x ^(r k , n ) :=J2d) {ei ~ e2 ~ K) < 1010|X0 > W 2 < £i0e 2 0\K0 > { ^ 
if 



m.n 

With this definition we write the matrix element as 

{{t I 1 



(3.11) 



(3.12) 



[(-) fe "^ ZE <^ll^ Al H^l>] [(-) fe *- ^^<^1||[^ 2 ^] A2 ||P>] [(-) fe - ^ ZE <f> 2 ||[^3^] A |l^2>] 
Ai A2 A 

Rh(fk)Rhi (rkjK^ 2 ^ 1 ^ 2 (r k , ri)R Pl (ri)Rp(ri)K x,e ' i ^ 3 (r k , fj)Rp 2 { r j)Rh 2 ( r j) 
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'til X 2 i}^i°^0|AO)(-) A2 +' 



[(-) fc3 ' 2 ^ z: (P2||^ A ||/ l2 >] [(-) fe "^ ZL <^ll[^ 2 ^] Al |l^i>] [(-) fe3 ' 1 ^ !E <Pill[^3^] A2 lb>] 

A Ai A2 

Bp s (r*)\(r^ Al ' <1 '^rk,ri)4 1 (ri)4(r»)^ !A ' i 1r fe ,r j )4(r j )^(r i ) 

A 4 1 \„.n,.nu.nw xAx+fe ( 3 ' 13 ) 



+ U Ai A2 }^iO^O|A 2 0)(-) A - 



A2 A Ai 

R Pl (r k )R p (r k )K x ^ ^ (r k ,n)R hl ( n )R h {n)K x ^ (r k , r 3 )R P2 (r 3 )R h2 (r 3 ) 

Similarly, we get the commutator contribution as 

v tni,cx,X _ Q a A iA 2 \(fc fe + fc B1 +fc hl )f Ai A 2 Al 





1 


A 2 ^ 


«! 


1 


A 




1 


Ai J 



[(-) fcA ^<'i||[^a-] A MI'»i>] [(-^^^illP^^l^llP)] [(-^"^(PallP^^Hfta)] 
Ai A 2 A 

r4 1 An 

^ L 1 A > (3.14) 

[(-D^T^IIP^o-^Hfea)] [(-)^^<ft||p^ CT ]^||/n>] [(-)^ 1 ^<Pi||[^.^] Aa lb>] 
A Ai A 2 

R P2 ( r k)Rh 2 ( r k)K Xueue2 (r k ,r i )R hl (r i )R h (r i )K X2 ' i4 ' i:i (r k ,r : j)R pi (r 3 )R p (r 3 ) 

(h 1 X) 

+ < h 1 Ai \\ 2 (£iO£±0\LO)(-) X2+l2+l3 
II 1 A 2 J 



[C-)**' i: ^ Z: <l»xllP^^] Aa IIJ»>] [(-)^= : ^ Z[: (Pa|IP^»^] A ll^>] [(-^^^^IIP^.^MI'H)] 
R Vl {r h )R p {r h )K XM {r k ,r i )R h M)Mn)K w ^ 
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